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ABSTRACT 

The 09IV star HD 57682, discovered to be magnetic within the context of the MiMeS 
survey in 2009, is one of only eight convincingly detected magnetic O-type stars. 
Among this select group, it stands out due to its sharp-lined photospheric spectrum. 
Since its discovery, the MiMeS Collaboration has continued to obtain spectroscopic and 
magnetic observations in order to refine our knowledge of its magnetic field strength 
and geometry, rotational period, and spectral properties and variability. In this paper 
we report new ESPaDOnS spectropolarimetric observations of HD 57682, which are 
combined with previously published ESPaDOnS data and archival Ha spectroscopy. 
This dataset is used to determine the rotational period (63.5708 ± 0.0057 d), refine 
the longitudinal magnetic field variation and magnetic geometry (dipole surface field 
strength of 880 ±50 G and magnetic obliquity of 79 ±4° as measured from the magnetic 
longitudinal field variations, assuming an inclination of 60°), and examine the phase 
variation of various lines. In particular, we demonstrate that the Ha equivalent width 
undergoes a double- wave variation during a single rotation of the star, consistent 
with the derived magnetic geometry. We group the variable lines into two classes: 
those that, like Ha, exhibit non-sinusoidal variability, often with multiple maxima 
during the rotation cycle, and those that vary essentially sinusoidally. Based on our 
modelling of the Ha emission, we show that the variability is consistent with emission 
being generated from an optically thick, flattened distribution of magnetically-confined 
plasma that is roughly distributed about the magnetic equator. Finally, we discuss our 
findings in the magnetospheric framework proposed in our earlier study. 

Key words: stars: individual HD 57682, stars: magnetic fields, stars: circumstellar 
matter, stars: rotation, stars: winds, outflows, techniques: polarimetric 



* Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research 
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1 INTRODUCTION 



HD 57682 is an 09 subgiant, and one of only eight 
O-type stars convincingly detected to host magnetic 
fields: the zero-age main sequence (ZAMS) 07V star 1 
Ori C=HD 37022 (jPonati et al.1 120021 : IWade et all l2006h 
the more evolve d Of?p stars HP 108 (iMartins et"all 
l2010h , HD 1 48937 faubrig et al.ll2008l: IWade et alJl2012aF 7 



HD 191 612 dDo"nati et al.1 120061 7 IWade et al.1 120 llh. NGC 
1624- 2 (IWade et al.l|2012bh and CPD -28 2561 (jHubrig et al.1 
1201 ll : IWade et al.1 l2012cf). the re c ently detected Carina 
nebula star Trl6-22 ( (Naze et al.1 l2012h . and HD 57682 
(jGrunhut et al.1 l2009h . In addition, a small number of 
other O-type stars have been tentatively reported to 
be mag netic in modern literature (e.g . iBouret et al.l 



2008: Hubrig et al 



Hubrig et all l201lh 



l2008l : lHubrig. Oskinova fe Schollerll2oTTI : 

These stars have either been found 



through indepe ndent observation and /or analysis to be non - 
magnetic (e.g. iFullerton et~ai1 l201ll : iBagnulo etld1l2012h . 
or have yet to be independently re-observed or re- analysed. 
These small numbers are both a reflection of the rarity of O- 
type stars with detectable magnetic fields, and the challenge 
of detecting such fields when present. 



Based on seven spec t ropo larimetric observations of 
HD 57682, iGrunhut etaD (|2009h characterised its physi- 
cal and wind properties, deriving, in particular, a mass 
of 17lg 9 M , a radius of 7.0tigi?©, a mass- loss rate 
logM — — 8.85M0yr -1 and wind terminal velocity Voo — 
1200^200 kms -1 . Their spectroscopic and magnetic data in- 
dicated that the star was variable, likely periodically on a 
time- scale of a few week s . Usi n g the Bayesian statis tical 
method of IPetit fe Wade] (|2012h . IGrunhut et al.1 (|2009h es- 
timated the surface dipole strength to be I68OI355 G. This 
magnetic field strength, combined with the inferred physi- 
cal and wind parameters, indic ated a wind magnetic con - 
finement parameter of 1.4 x 10 4 ([ud-Doula &; Owockill200"2h , 
leading the authors to suggest that the observed line profile 
variations were a consequence of dense wind plasma confined 
in closed magnetic loops above the stellar surface. 
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of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii, 
f E-mail: Jason. Grunhut @rmc.ca 



The present study seek s to refine and elabora te the pre- 
liminary results reported bv I Grunhut et al . (2009). In Sect. [2] 
we discuss the new and archival observations used in our 
analysis, including extraction of Least-Squares Deconvolved 
(LSD) mean profiles and measurement of the longitudinal 
magnetic field. In Sect. [3] we perform a period analysis of 
the Ha and longitudinal field data and demonstrate that a 
single period of ~63 d is required to reproduce the periodic 
variations of the measurements, which span over 15 years. 
In Sect. [4] we discuss the implications of the inferred rota- 
tion period for the projected rotational velocity of the star, 
and, in particular, th e estimate of 15 kms -1 proposed by 
IGrunhut et al.l (|2009h . In Sect. [5] we constrain the magnetic 
field properties using the longitudinal field measurements 
and direct fits to the observed LSD Stokes V profiles. In 
Sect. [6] we present an extensive overview of the variability 
of line profiles of H, He and metals. In Se ct. |7] we revisit 
the m agnetospheric framework proposed bv IGrunhut et all 
(2009), and in Sect. [8] we discuss our results and present our 
conclusions. 
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2 OBSERVATIONS 

2.1 Polarimetric data and measurements 

A total of 38 polarimetric observations were collected with 
the high-resolution (R ~ 68 000) ESPaDOnS spectropo- 
larimeter mounted on the Canada-France-Hawaii telescope 
(CFHT). The observations were obtained in the context of 
the Magnetism in Massive Stars (MiMeS) CFHT Large Pro- 
gram on 20 different nights over a period of two years from 
2008 December to 2010 December. Each observation consists 
of a sequence of four sub-exposures that were processed us- 
ing th^_UPENAnipej^ne running libre-esprit, as described 
bv lDonati et al.l (|1997). The final reduced products consist 
of an unpolarised Stokes I and circularly polarised Stokes V 
spectrum. The individual sub-exposures are also combined 
in such a way that the polarisation should cancel out, re- 
sulting in a diagnostic null spectrum that is used to test 
for spurious signals in the data. During most nights two po- 
larimetric sequences were obtained and the un-normalised 
spectra were combined to produce 20 distinct observations. 
Se ven of our twenty ob servations were previously discussed 
bv lGrunhut et al.l (|2009h and we present 13 new observations 
here, the details of which are listed in Table [U In Fig. [T] we 
show spectra for selected spectral regions at a few differ- 
ent rotational phases to highlight the variability observed in 
almost all spectral lines (see Sect. [3] for further details). 

From each observation we extracted LSD pr ofiles using 
the technique as described bv lDonati et al.l (|l997l ). The mea- 
su rements reported here differ in details from those reported 
bv lGrunhut et al.1 12009) as we have adopted an updated line 
mask with significantly more lines, but the measurements 
are fully consistent within their respective uncertainties. Our 
new mask has an average Lande factor and average wave- 
length of 1.180 and 4682 A versus 1.135 and 4649 A from the 
original mask. The new mask w as derived from the V ienna 
Atomic Line Database (VALD; Pis kunov et al.|[l995h for a 
^eff = 35 000 K, log(g) = 4.0 model atmosphere. Only lines 
with intrinsic depths greater than 1 percent of the contin- 
uum were included in the line list over the full ESPaDOnS 
spectral range (370 to 1050 nm), resulting in 1648 He and 
metallic lines. The line list was further reduced by remov- 
ing all intrinsically broad He and H lines, all lines that are 
blended with these lines, in addition to all lines blended with 
atmospheric telluric lines. This resulted in a final list of 571 
metallic lines used to extract the final LSD profiles com- 
puted on a 1.8kms _1 velocity grid, as illustrated in Fig. [2] 
The detection probability for eac h LSD spectrum wa s com- 
puted according to the criteria of iDonati et al.l (|l997h . Each 
observation resulted in a definite detection (DD) of an ex- 
cess signal in Stokes V within the line profile (False Alarm 
Probability (FAP) < 10 -5 ), except for the observations ob- 
tained on the nights of 2008-12-05 and 2009-05-07, which 
were non-detections (ND). 

We calculated the longitudinal magnetic field (Be) 
and null measurements (Ne) from each LSD profile using 
the fir st-order moment method discussed bv lRees &; Semell 
(| 19791 ). using an integration range between -25 and 
75kms _1 (see Tabled]). This range was adopted to encom- 
pass the full profile range of the Stokes V signature. The 
resulting Be measurements range from -170 to -350 G with 
a median uncertainty of 20 G. 



2.2 Spectroscopic data and measurements 

In addition to the polarimetric spectra, we also utilised 
archival spectroscopic observations to analyse the line pro- 
file variability (LPV) of Ha. Seven spectra were obtained 
with the Coude Auxiliary Telescope (CAT) using the Coude 
Echelle Spectrograph (CES, R ~ 45 000) at the European 
Southern Observatory (ESO), La Silla, Chile in 1996 Febru- 
ary. These observations were obtained as part of a large 
moni toring campaign of bright O-type stars, as described 
by iKaper et al.l (Jl998). Eight spectra were obtained with 
the Ultraviolet and Visual echelle Spectrograph (UVES, 
R rsj 50,000) mounted on ESO's Very Large Telescope 
(VLT) on 2002-12-13. These spectra were combined to yield 
a single observation. Five spectra were obtained with the 
Fibre- fed Extended Range Optical Spectrograph (FEROS, 
R rsj 48, 000) mounted on ESO's 2.2-m telescope at La Silla: 
four observations between 2004-02-07 and 2004-02-10 (on 
three separate nights), and another observation on 2004- 
12-22. Between 2009-10-07 and 2009-10-09 ten spectra were 
obtained with the Echelle spectrograph (R ~ 40, 000) on 
the 2.5-m Du Pont telescope at the Las Campanas Observa- 
tory (LCO). The nightly spectra were combined to produce 
three distinct observations. Lastly, we also used spectra from 
the BeSS database obtained between 2008 November and 
2011 April. Fifty-one spectra were collected, primarily with 
a LHIRES3 spectrograph (R ~ 15 000) using various tele- 
scopes with diameters between 0.2 to 0.3-m. A log of these 
spectroscopic observations is listed in Table [2] 

From the spectroscopic datasets we measured the equiv- 
alent width (EW) variations of Ha (in addition to several 
additional spectral lines present in the ESPaDOnS dataset - 
see Sect. [6] for f urther details), follow ing the same procedure 
as described bv I Wade et all (|2012aT ). The Ha EW measure- 
ments are listed in Table [3] 

We also utilized archival IUE spectra, the details of 
which are described bv lGrunhut et al.1 (|2009h . 



2.3 Photometric data 

Photometric observations were also utilised to study the op- 
tical variability (see Sect. [3]) and to characterise the stellar 
and circumstellar emission properties (see Sect. 17.2)) . 

The optical photometric ti me series corresponds to 85 
measurements from Hipparcos (Per rvman et alJ ll997) col- 
lected on 64 different nights between 1990-04-03 and 1993- 
06-06. 

Additional photometric observations used to construct 
the spectral energy distribution (SED) consists of optical 
data obtained from the Centre de Donnees astronomiques 
de Strasbourg (CDS), and infrared measurements from the 
Two Micron All Sky Survey (2MASS), the Deep Near In- 
frared Survey (DENIS), the Infrared Astronomical Satel- 
lite (IRAS), the AKARI satellite and the Wide-field In- 
frared Survey Explorer (WISE), which were all obtained 
from NASA's Infrared Science Ar chive (IRS A). Additionally, 
we also used photometry from iThe. Wesselius &; Janssenl 
(| 19861 ). 

The optical and infrared magnitudes were converted to 
flux measurements at the zero-point of the associated pho- 
tometric filters using conventional methods. 
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Table 1. Journal of polarimetric observations listing the date, the heliocentric Julian date (2 450 000+) , the number of sub-exposures 
and the exposure time per individual sub-exposure, the phase according to Eq. [TJ the peak signal-to-noise (S/N) ratio (per 1.8 kms -1 
velocity bin) in the V-band , the m ean S/N in LSD Stokes V profile, the evaluation of the detection level of the Stokes V Zeeman signature 
according to iDonati et al.l (|l997T ) (DD=defmite detection, ND=no detection), and the derived lon gitudinal field and lon gitudinal field 
detection significance z from both V and N. The first seven observations were already discussed by I Grunhut et all (2009). 





HJD 


^exp 




PK 


LSD 


Det. 


V 




N 




Date 


(2 450 000+) 


(s) 


Phase 


S/N 


S/N 


Flag 


B £ ± * B (G) 


z 


N £ ± a N (G) 


z 


2008-12-05 


4806.0798 


1 x 4 x 500 


0.9408 


311 


2186 


ND 


354 ± 93 


3.8 


-59 ± 93 


0.6 


2008-12-06 


4807.1081 


1 x 4 x 500 


0.9570 


949 


8339 


DD 


248 ± 24 


10.3 


-11 + 24 


0.5 


2009-05-04 


4955.7675 


1 x 4 x 540 


0.2955 


1463 


12769 


DD 


-42 ± 17 


2.5 


-19 + 17 


1.1 


2009-05-05 


4956.7498 


2 x 4 x 600 


0.3109 


1388 


12192 


DD 


-64 ± 18 


3.6 


-23 + 18 


1.3 


2009-05-07 


4958.7805 


2 x 4 x 540 


0.3429 


653 


5546 


ND 


-78 ± 39 


2.0 


-17 + 39 


0.4 


2009-05-08 


4959.7489 


1 x 4 x 540 


0.3581 


1094 


9515 


DD 


-153 + 23 


6.7 


9 + 23 


0.4 


2009-05-09 


4960.7480 


2 x 4 x 540 


0.3738 


900 


7634 


DD 


-126 + 29 


4.4 


-7 + 29 


0.2 


2009-12-31 


5197.0957 


2 x 4 x 540 


0.0917 


1566 


14309 


DD 


233 ± 14 


16.2 


-7 + 14 


0.5 


2010-01-03 


5201.0427 


2 x 4 x 540 


0.1538 


1110 


10354 


DD 


133 ± 20 


6.6 


-1 + 20 


0.0 


2010-01-23 


5219.9545 


2 x 4 x 540 


0.4512 


1376 


12569 


DD 


-169 + 18 


9.6 


-13 + 18 


0.8 


2010-01-29 


5225.9836 


2 x 4 x 540 


0.5461 


1610 


14758 


DD 


-174 + 15 


11.5 


12 + 15 


0.8 


2010-01-31 


5228.0015 


2 x 4 x 540 


0.5779 


1179 


10817 


DD 


-172 + 20 


8.5 


-33 ± 20 


1.6 


2010-02-01 


5228.9032 


2 x 4 x 540 


0.5920 


1250 


8088 


DD 


-162 + 27 


6.0 


-8 + 27 


0.3 


2010-02-24 


5251.9131 


2 x 4 x 414 


0.9540 


1198 


10940 


DD 


234 ± 18 


12.7 


34 + 18 


1.9 


2010-02-28 


5255.9891 


2 x 4 x 414 


0.0181 


932 


7208 


DD 


190 ± 28 


6.8 


-33 ± 28 


1.2 


2010-03-04 


5259.9420 


2 x 4 x 414 


0.0803 


1163 


10477 


DD 


208 ± 20 


10.6 


-10 + 20 


0.5 


2010-03-08 


5263.9096 


2 x 4 x 414 


0.1427 


878 


7921 


DD 


187 ± 26 


7.1 


-4 + 26 


0.2 


2010-11-28 


5529.0352 


2 x 4 x 415 


0.3133 


1224 


10980 


DD 


-60 + 19 


3.1 


-5 + 20 


0.3 


2010-12-24 


5555.0639 


2 x 4 x 415 


0.7227 


1083 


9567 


DD 


-13 + 22 


0.6 


-27 + 22 


1.2 


2010-12-30 


5561.0731 


2 x 4 x 415 


0.8172 


995 


8850 


DD 


89 + 23 


3.8 


20 + 23 


0.9 



Table 2. Journal of spectroscopic observations listing the iden- 
tification of the dataset, the typical resolving power of the in- 
strument, the epoch of the observations, the number of spectra 
obtained within the given dataset, the resulting number of dis- 
tinct observations, and the median per pixel signal-to-noise (S/N) 
ratio in the Ha region for the dataset. 





Res. 




Num 


Num 


Median 


Name 


Power 


Epoch 


Spec 


Obs 


S/N 


CES 


45 000 


1996 


7 


7 


540 


UVES 


50 000 


2002 


8 


1 


160 


FEROS 


48 000 


2004 


5 


5 


170 


LCO 


40 000 


2009 


10 


3 


360 


BeSS 


15 000 


2008-2011 


51 


51 


70 



Table 3. Spectroscopic Ha EW measurements. Included is the 
instrument or observatory name, the heliocentric Julian date, and 
the measured Ha equivalent width (EW) and its corresponding 
la uncertainty. Shown here is only a sample of the table. The full 
table can be found in Appendix [Al 



Dataset HJD Ha EW (A) a E w (A) 



CES 


2450125.6216 


1.270 


0.014 


CES 


2450126.6291 


1.108 


0.007 


CES 


2450127.6296 


0.926 


0.008 


CES 


2450128.6284 


0.735 


0.012 


CES 


2450129.6287 


0.569 


0.011 


CES 


2450130.6194 


0.454 


0.010 



3 PERIOD ANALYSIS 

I Grunhut et al . (2009) attempted to infer the rotational pe- 
riod (Prot) of HD 57682 based on the limited temporal sam- 
pling of their dataset. With the larger dataset presented in 
this paper, we re-evaluate the periodicity and rotation pe- 
riod of HD 57682. 

Using a Lomb-Scargle technique ([Press et al.lll996l ). we 
analysed the longitudinal magnetic field variations Bg and 
the EW measurements of Ha. The periodograms obtained 
from these analyses are displayed in Fig. [3] The results ob- 
tained from the Bg measurements show a clear peak in the 
periodogram at a period of 63.521°;^ d, where the uncer- 
tainties represent the la limits. When phased with this pe- 
riod, the Bi measurements show a clear sinusoidal varia- 
tion. However, the periodogram obtained from the Ha mea- 
surements shows a clear peak occurring at 31.7929io'.ooi4 °^ 
roughly half the period that is obtained from the Bt mea- 
surements. If we phase the Ka EW measurements to this 
period, we obtain a roughly sinusoidal variation, but with a 
relatively broad dispersion. On the other hand, the Bg mea- 
surements cannot be reasonably phased with this period. We 
proceeded to obtain a periodogram o f the Ha data using the 
multi- harmonic fitting technique of Sch warzenberg-Czernvl 
( 1996), including contributions from the first harmonic only. 
The resulting periodogram shows two peaks, one occurring 
at about 31 d and the other occurring at 63.5708 =b 0.0057 d. 
This second peak is consistent with the period found from 
the magnetic measurements. Therefore, adopting this period 
as the stellar rotation period and the phase of maximum 
positive Bg as HJDo we derive the following ephemeris: 

HJD^ = 2453347.71(35) + 63.5708(57) • E, (1) 
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Figure 1. Selected regions of visible spectra from ESPaDOnS. Spectra at different rotational phases are shown to highlight the observed 
variability in almost all lines (0.02 - solid black, 0.30 - dashed red, 0.45 - dotted blue, 0.72 - dash-dotted green; see Sect. [6] for further 
details). The ions with the most significant contribution to each line are also labelled. 



where the uncertainties (la limits) in the last digits are in- 
dicated in brackets. Unless otherwise stated, all further data 
are phased according to this ephemeris. 

In Fig. [4] we show the phased Bg measurements, Ha 
EW variations, and Hipparcos photometry of HD 57682. The 
Bi measurements follow a simple, sinusoidal variation with 
maximum positive Bg occurring at phase 0.0 and maximum 
negative Bi occurring at phase 0.5. The Ha EW measure- 
ments show a double peaked variation with maximum emis- 
sion occurring at phases 0.0, and 0.5, and minimum emission 
occurring at phases 0.25 and 0.75. This is further discussed 
in Sect. [7] The Hipparcos measurements do not show any 
significant variations when phased with the rotation period. 
Furthermore, the periodogram shows no significant power at 
this or any other period and the scatter of the photometric 
variations and their associated uncertainty is fully consistent 
with no variability. 



4 ROTATIONAL BROADENING 

iGrunhut et al.l (|2009h measured the projected rotation al ve- 
locity (v sin i) using the Fourier transform metho d (e.g. lGravl 
ll98ll : IJankoJl995inSim6n-Dfaz fe Herrerdl2007h . They com- 
pared the positions of the first nodes in the Fourier spectrum 
to a theoretical Gaussian profile convolved with a rotation- 
ally broadened profile corresponding to a particular vsini, 
and found a vsini = 15 =b 3kms -1 . 

If we assume rigid rotation and take the radius as in- 
ferred by IGrunhut et all (|2009h (R* = 7.0±l£ Rq) and the 
period obtained in this work (P = 63.5708 d), we find a max- 
imum allowed vsini (i.e. for sini = 1) ^7.5 km s -1 , which 
is inconsistent with the results obtained from the Fourier 
analysis. This implies th at contributions from n on- rotational 
broadening, as noted by IGrunhut et al.l {2009), significantly 
affect the Fourier spectrum resulting in the incorrect inter- 
pretation of the first node representing the rotational broad- 
ening. 
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Figure 2. Mean LSD Stokes V (top), diagnostic null (middle) 
and Stokes / profiles (bottom) of HD 57682 from observations ob- 
tained at four different rotational phases (see Sect. [3] 0.02 - solid 
black, 0.30 - dashed red, 0.45 - dotted blue, 0.72 - dash-dotted 
green). The V and N profiles are expanded by the indicated fac- 
tor and shifted upwards, and smoothed with 3 pixel boxcar for 
display purposes. A clear Zeeman signature is detected at most 
phases in the Stokes V profiles, while the null profile shows no 
signal. The integration limits used to measure the longitudinal 
field are indicated by the dotted lines and the mean la pixel 
uncertainty for each profile is also shown. 



o 

Oh 



P=31.7929d 



0.02 0.03 

Frequency (1/d) 

Figure 3. Periodograms obtained from the Bg measurements 
(solid black) and from the Ha EW variations including contri- 
butions from the first harmonic (dashed red). Note the strong 
power present in the Ha periodogram at P = 31.7927 d but the 
lack of power from the B% measurements at this period. However, 
a consistent period is found at P = 63.5708 d from both datasets. 



In our endeavour to constrain the rotational broaden- 
ing, we re-analysed the spectra, particularly focusing on 
those obtained at phases where the lines appeared their nar- 
rowest and therefore were least affected by macroturbulence 
(see Sect. [6] for furt her details). We identified lines from a 
theoretical synth3 (jKochukhovl l2QQ7h spectrum for which 
the resulting Fourier spectrum was minimally affected when 
significant contributions of macroturbulence were included. 
We then compared the Fourier spectrum obtained from a 
theoretical profile with 5kms _1 rotational broadening to 



that obtained from the observed spectrum and measured 
the difference between the positions of the first nodes to ob- 
tain the vsini of the spectral line. Using this procedure on 
several lines from the spectrum obtained on 2009 December 
31 indicate a vsini — 6.1 =b 1.9kms _1 , consistent with the 
range implied by the rotation period and estimated radius of 
HD 57682. However, we note that this period is sufficiently 
imprecise that it cannot be used to usefully constrain the 
inclination of the rotation axis. We also found that the C IV 
A5801 line was one of the lines least affected by the addition 
of macroturbulence and therefore also used this single line 
to determine a unique vsini, as illustrated in Fig. [5] From 
measurements of this line at several rotation phases we find 
a best-fit vsini — 4.6 =b 0.6kms _1 , where the uncertainties 
represent the la limits. Given the range of vsini obtained 
from the C IV A5801, our adopted R* and P ro t, implies that 
i rsj 56^23° 5 or that we can place a lower limit on i such 
that i > 30°, at a la confidence. However, at ~4.6kms -1 
we are approaching the spectral resolution of ESPaDOnS 
(Mikrns -1 ) and therefore one must be cautious to not 
over interpret these results. 



5 MAGNETIC FIELD GEOMETRY 

The magnetic field geometry of HD 57682 was investigated 
assuming the field is well described by the Oblique Rotator 
Model (ORM), which is characterised by four parameters: 
the phase of closest approach of the magnetic pole to the 
line-of-sight 0o, the inclination of the stellar rotation axis i, 
the obliquity angle between the magnetic axis and the ro- 
tation axis P and the dipole polar strength Bd. This model 
naturally predicts the simple sinusoidal variation of the lon- 
gitudinal magnetic field measurements as shown in Fig. 2] 
implying that the longitudinal field curve of HD 57682 is well 
described by a simple dipole topology of its global magnetic 
field. 

We first attempted to infer the characteristics of the 
magnetic field based on the Bg measurements. These mea- 
surements are nearly symmetric about and we can there- 
fore infer that either i or ft (or both) are clos e to 90° . This 
is sim ilar to the Bg variation of f3 Cephei (e.g. iDonati et al.l 
2001), and similar to this star, Bd and f3 are well-constrained 
except as i approaches or 90 °. As i approaches 0°, Bd be- 
comes unconstrained, while as i approaches 90°, /3 becomes 
unconstrained. To infer the characteristics of the magnetic 
field, we carried out a \ 2 minimisation, comparing the ob- 
served Bi curve to a grid of computed longitudinal field 
curves to determine Bd and (3 for a fixed linear limb dark- 
ening coefficient of 0.35 (|Claretll2000h and various fixed in- 
clinations. 

The resulting \ 2 distributions from our fits are shown 
in Fig. [6] For an assumed i = 60° we find that Bd — 880 =b 
50 G and f3 = 79 ± 4°. If we take i = 30°, we infer that 
Bd = 1500 =b 90 G and f3 = 86 =b 2°, which places an upper 
limit on the dipole field strength if i > 30° as implied from 
the constraints derived in Sect. [4] Also shown in Fig. [6] are 
X 2 distribution for other inclinations to illustrate how the 
distribution varies as a function of inclination. 

We also attempted to constrain the magnetic field char- 
acteristics based on fits to the individual LSD Stokes V pro- 
files by comparing our profiles to a large grid of synthetic 
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Figure 4. Phased observational data using the ephemeris of Eq. ^ Upper: Longitudinal magnetic field variations measured from LSD 
profiles of the ESPaDOnS spectra. Middle: Ha equivalent width variations measured from ESPaDOnS (black circles), CES (red squares), 
FEROS (blue Xs), Las Campanas Observatory (turquoise diamonds), UVES (pink up-facing triangles), and BeSS (green down-facing 
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data. 
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Figure 5. Fourier spectrum of the C IV A5801 line (black) com- 
pared to a theoretical line profile with vsini = 5kms _1 with 
an additional macroturbulent velocity of 20kms _1 (red), which 
we find necessary to provide a good fit to the observed line pro- 
file. The difference in the positions of the first nodes (as indi- 
cated by the vertical lines) are used to determine the best-fit 
vsini = 4.2kms — 1 for this observation. 



profiles, also characterised by the ORM. The models were 
computed by performing a disc integration of local Stokes V 
profiles assuming the weak field approximation and uniform 
surface abundance. The parameters of the Stokes V profiles 
were chosen to provide the best fit to the width of the av- 
erage of all the observed / profiles, and the best fit to the 
depths of each individual observed I profile. For each ob- 
servation we found the parameters that provided the lowest 

© 2012 RAS, MNRAS 00Q.[TU2T1 



X 2 value and inferred the maximum likelihood by combining 
these results in a Bayesian framework. 

We conducted our fits for varying inclinations and found 
that the maximum likelihood model for i ~ 80° provided 
the lowest combined total % 2 5 but this value was not signif- 
icantly different from those obtained for other inclinations. 
However, as i increased we did find increasing deviations 
between the best 2^-fit and profile-fit model parameters; at 
i = 30° we found a best i^-fit f3 = 86° and a profile-fit 
P = 80°; at % = 80° we found a best B e -fit f3 = 56° but 
P — 34° from the profile fits. We therefore conclude that we 
are unable to constrain the inclination from fits to the LSD 
profiles and continue to adopt i — 60° for the remainder of 
this discussion as this is the inclination inferred from the 
vsini measurements and provides similar magnetic param- 
eters from both the Bi and profile fits. 

The maximum likelihood profile-fit model for i — 60° 
was found with Ba = 700 =b 30 G, where the uncertainty 
corresponds to the 95.6 percent significance level, and f3 — 
68 =b 2°, with the uncertainty corresponding to the 98.2 per- 
cent confidence limits. While the formal uncertainties are 
very low, we note that individual best-fit parameters to each 
observation range from B d - 400-2000 G and ft ~ 0-150°. 
However, while the fits to Bd and f3 are not well constrained 
from fits to the individual observations at most phases, the 
maximum likelihood model is well constrained. 

In Fig. [7] we compare the observed Stokes V profiles 
(grey circles) with the best profile-fit model for the given ob- 
servation (dash-dotted blue), the maximum likelihood LSD 
model (dotted red) and the model corresponding to the best- 
fit parameters that were obtained from the Bg measurements 
for % = 60° (B d = 880 G and f3 = 79°; solid green). The 
quality of the fits of the best-fit and maximum likelihood 
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Figure 6. Reduced x 2 contours of dipole field strength versus magnetic obliquity (5 permitted by the longitudinal field variation of 
HD 57682 for varying inclinations. 



models are similar, which shows that a single dipole con- 
figuration can reasonably reproduce the observed Stokes V 
profiles. However, systematic differences, particularly in the 
wings of the profiles, are evident between the models and the 
observations, resulting from the fact that our simple model 
cannot reproduce the extended wings in the observed I or 
V profiles. It is also evident that the model corresponding 
to the best-fit to the Bg curve does a poorer job of fitting 
the observed Stokes V profiles at several phases, possibly 
indicating that there may be an issue with the way we chose 
to model the line profiles. 



We investigated this potential issue by re-conducting 
our fits, but neglecting the varying line depths for each ob- 
servation as was done in the original procedure. Instead we 
adopted a line depth as determined from the average of all 
our observations. From the new fits for i = 60° we found 
a best-fit Bd ~ 900 G, which is consistent with the results 
of the Bg measurements. However, we find a larger disagree- 
ment between the inferred magnetic obliquity as f3 was found 
to be ^60°. We also note that the overall combined total % 2 
value was slightly worse from these fits. 



Another possibility is that the global magnetic topology 
contains significant contributions from higher order multi- 
poles. The Bg curve is rather insensitive to these higher order 
moments, in contrast to the velocity-resolved Stokes V pro- 
files. We investigated this possibility by applying our fitting 
procedure to a timeseries of synthetic profiles calculated as- 
suming our original dipole, supplemented by a significant 
aligned quadrupolar field moment. This model yields a lon- 
gitudinal field curve essentially indistinguishable from that 
produced by the pure dipole. Our tests indicate that by try- 
ing to fit a pure dipole model to the dipole+quadrupole pro- 
files, the procedure would systematically infer dipole field 
strengths that were strongly incompatible with the results 
obtained from the Bg curve. Specifically, our tests showed 
a systematic over- estimation of the dipole field strength at 
most phases. This trend is inconsistent with the measure- 
ments presented in Fig. [71 and we therefore conclude that 
our polarimetric observations do not suggest any significant 
detectable contributions from higher order multipoles. 



6 LINE VARIABILITY 

6.1 Equivalent width and radial velocity 
measurements 

HD 57682 stands out from the sample of magnetic O-type 
stars due to its sharp spectral lines and the apparently lim- 
ited contamination of its spectrum by wind and circumstel- 
lar plasma emission. D etection and measu rement of the line 
profile variations (e.g. IStahl et "all [l996) can therefore be 
performed more sensitively for a much larger sample of lines 
than for any other magnetic O-type star. It therefore repre- 
sents a uniquely suited target for understanding the general 
influence of magnetised winds of O-type stars on their spec- 
trum and spectral variability. 

To quantify the variability, we used the EW variations 
for numerous spectral lines. We also measured the radial ve- 
locities of each spectral line by fitting a Gaussian to the cores 
of the profiles to determine the central wavelength at each 
phase. A number of other methods were investigated, but we 
ultimately settled on fitting the cores using a Gaussian as 
it provided less scatter in the measurements. However, since 
Ha displayed too much emission to carry out this fitting, we 
quantified the radial velocity of Ha by measuring the centre 
of gravity (or centroid) of the line profile. 

Our observations reveal that HD 57682 displays a rich 
variety of line profile variability. In addition to the double- 
wave variations of the Ha line discussed in Sect. [3] we ob- 
serve variations in essentially every line present in the stellar 
spectrum, including other H Balmer lines, lines of neutral 
and ionised He, and lines of metals: Cm and iv, Niii, On, 
Mgn, and Sim and IV. The EW and radial velocities of all 
lines vary coherently when phased according to the rota- 
tional ephemeris (Eq. [T]). 

The line variations can be approximately grouped into 
two categories. The first represents lines (mostly metals and 
some neutral He lines) with EW measurements that vary 
sinusoidally, as displayed in Fig. [8] The second includes 
those, like Ha, that display EW phase variations that are 
distinctly non-sinusoidal (as shown in Fig. [9}, and likely di- 
rectly reflect, to varying degrees, the variable projection of 
the flattened equatorial magnetospheric plasma density en- 
hancement, as further discussed in Sect. [7] These groups are 
discussed in more detail below. 

Except for Civ A5801,5812, all lines with sinusoidally 
varying EW measurements show minimum emission (or 
maximum absorption) at about phase 0.0 and maximum 
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Figure 7. Mean circularly polarised LSD Stokes V profiles (grey circles). The error bars represent the la uncertainties of each pixel. 
Also shown are the individual best-fit model profiles for each phase (dash-dotted blue), the model that provides the global maximum 
likelihood (dotted red; = 700 G, (3 = 68°) and the model corresponding to the best-fit to the Bg measurements (B d = 880 G, (3 = 79°; 
solid green), all for an assumed inclination of 60°. The phase and best-fit parameters of the individual observations are also indicated. 



emission (or minimum absorption) at phase 0.5. At phase 
0.5, these lines are shallower and broader than at phase 0.0, 
as illustrated in Fig. [T] The difference in line width and 
depth is quite significant. Most lines have an extended red 
wing at phase 0.5 that disappears at phase 0.0. 

The Ha line shows variable emission with two max- 
ima per rotation cycle. The phases of the emission extrema, 
0.0 and 0.5, correspond to the phases of longitudinal mag- 
netic field extrema and therefore closest approach of the 
magnetic poles to the line- of- sight. In addition, the Ha EW 
curve exhibits two emission minima at phases 0.25 and 0.75, 
corresponding to the phases where the magnetic equator 
crosses the lin e- of- sight. These p henomena are usually in- 
terpreted fe.g. lDonati et al. 2001) as the consequence of the 
variable projection of a flattened distribution of magneto- 
spheric plasma trapped in closed loops near the magnetic 



equatorial plane, and possibly the occultation of the stel- 
lar disc by this plasma at phases 0.25 and 0.75. Only three 
other lines - Hei A5876, HeiiA4686 and R/3 (see Fig. [9) - 
show clear evidence of similar reversals at phases 0.25 and 
0.75, although a marginal contribution is suspected for He I 
A4921 and Hei A4471. The emission level in Ha at phase 
0.5 is roughly 10 percent larger than at phase 0.0. Likewise, 
the emission minimum at phase 0.75 appears 20 percent less 
than the corresponding emission minimum at phase 0.25. 
However, these measurements are primarily measured from 
BeSS spectra with limited S/N. 

In almost all the spectral lines that show sinusoidally 
varying phased EW measurements, we find approximately 
sinusoidally varying radial velocity measurements. These 
variations reach a maximum radial velocity between phases 
0.2-0.3 and show a peak-to-peak amplitude of about 
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2kms _1 , with a mean radial velocity of ^22kms -1 . This 
is not the case for H7, which reaches a maximum radial ve- 
locity at phase 0.0 and minimum at phase 0.5, and shows 
a peak-to-peak difference of ~8kms _1 . Other than H/3, we 
see non-sinusoidal variations in the radial velocity measure- 
ments for the lines that show clear, non-sinusoidal EW varia- 
tions, which often mirror the EW variations with maximum 
radial velocity occurring at minimum emission and mini- 
mum radial velocity occurring at maximum emission. The 
radial velocity measurements for H/3 phase in a similar man- 
ner with H7 - they are nearly sinusoidal, reach maximum 
radial velocity at phase 0.0 and appear flat-topped. 

6.2 Line profile variations 

In Figs. [10] and [11] we plot the phased line profile varia- 
tions, or dynamic spectra. In these figures, we have sub- 
tracted the profile obtained on 24 December 2010 to high- 
light the variability. This spectrum was obtained at phase 
0.72 and corresponds to the observation where the magnetic 
equator is passing our line-of-sight and, as previously dis- 
cussed, when the flattened distribution of magnetospheric 
plasma is viewed nearly edge-on, and therefore provides the 
lowest contribution of emission. Our dynamic spectra indi- 
cate that a single pseudo-emission feature (a feature that 
appears in emission relative to the 24 December 2010 pro- 
file) occurs once per cycle, and is centred at phase 0.5 in 
the dynamic spectra of the lines with sinusoidally varying 
EW measurements (Fig. I10p . In almost all these lines, we 
find that the pseudo-emission is mostly confined to the in- 
ner core region and extends further out in the blue wing 
than the red wing. During phases of maximum emission we 
also find increased pseudo-absorption in the red wing, except 
for He 1 A4471, which does not appear to show additional 
pseudo-absorption, but does show strong evidence of pseudo- 
emission during phases of maximum absorption, which is 
centred at phase 0.0. Other lines show varying amounts of 
pseudo-emission during this phase as well, but this feature 
is most prominent in the He lines shown in Fig. 1101 

The dynamic spectra for the C IV lines appear different. 
Maximum emission still occurs once per cycle, centred at 
phase 0.5, but the emission feature is mainly constrained to 
the inner core region and does not extend into the blue wing 
as it does for other lines. During phases of maximum emis- 
sion, there appears to be an increase in pseudo-absorption 
in both wings, but the absorption in the red wing appears 
stronger. A similar dynamic spectrum is also found for the 
Hen A4686 line, as shown in the left panel of Fig. 1111 In 
almost all the lines shown in Figs. [10] and [TT] we find that 
the features all appear at low velocities with the centres 
of the features increasing in velocity with increasing phase. 
However, this does not appear to be the case with the He 
lines, which appear more symmetric, except for the extended 
pseudo-absorption into the blue wing, which is likely caused 
by the forbidden transitions in these He lines. 

In Fig. [12] we show the dynamic spectra of the other 
lines wit h non-sinusoidally vary ing EW variations. A NLTE 
tlusty f|Lanz & Hubenvl [2003) synthetic profile was sub- 
tracted from the Balmer lines in this figure to highlight the 
circumstellar emission contribution. As was found in the EW 
variations, we see two emission features per rotation cycle, 
one centred around phase 0.0 and the other centred around 



phase 0.5. In the Hei A5876, H7 and H/3 line we find the 
lines to be mainly in absorption and the emission feature at 
phase 0.0 to be narrower and weaker in emission than the 
feature at phase 0.5. In the He I line shown, both emission 
features appear redshifted with respect to the systemic ve- 
locity of HD 57682; the narrower feature is found to have a 
central velocity of ^75kms _1 while the broader feature has 
a central velocity of ^65kms -1 . There also appears to be 
an additional weak emission feature that is blueshifted with 
respect to the systemic velocity occurring at phase 0.0. H/3 
and H7 appear similar to the He I A5876 line. Both show 
a strong and broad emission feature at phase 0.5, and a 
weaker, narrower emission feature at phase 0.0. This is not 
the case for Ha, which shows both emission features with 
nearly the same width and same emission level; it is evident 
that the strength of the weaker emission feature at phase 
0.0 is decreasing for higher Balmer lines. Unlike the Hei 
A5876 line, the emission features of the Balmer lines have 
their central velocities alternating about the systemic veloc- 
ity; the feature at phase 0.0 appears blueshifted, while the 
feature at phase 0.5 is redshifted. 



7 CIRCUMSTELLAR ENVIRONMENT 
7.1 Magnetosphere 

As discussed in the previous section, the Ha emission vari- 
ation is often qualitatively attributed to the variable pro- 
jection of a flattened distribution of magnetospheric plasma 
trapped in closed loops near the magnetic equatorial plane. 
To this end, we explored the pote ntial of using a "Toy " 
model, similar to that described by iHowarth et all ([20071 ), 
to fit the observed Ha EW variations. This model assumes 
that the Ha emission is formed in a centred, tilted, infinitely 
thin, optically thick disc, in which the relative emission is 
only a function of the projected area of the disc, taking into 
account occult at ion by the star. We are only able to model 
the relative emission, but our model indicates that if we as- 
sume the inclination angle between the disc axis and the 
rotation axis (a) is the same as the magnetic obliquity (/?) 
then only subtle variations are predicted in the EW curve for 
different inclinations (assuming the relationship between i 
and P is constrained by the ORM), which are indistinguish- 
able with the precision of our current measurements. We 
also find that higher a values provide much better fits and 
therefore we adopt the a = j3 = 79° as derived from the Bg 
measurements as our disc inclination. Using this value, we 
first attempted to constrain the disc radii but realised that 
the inner (Ri n ) and outer (.Rout) radii of the disc are poorly 
constrained by the precision of our observations. However, 
a best-fit Ri n = 1.3 =b 0.3 i2* and R out — 1.8^0.2 ^* ( as mea_ 
sured from the centre of the star, i.e. a distance of 0.3±0.3 i?* 
and 0.8^0 2 from the stellar surface) are found, where the 
uncertainties represent the la limits. An illustration of this 
model is presented in Fig. 1131 

A comparison between the predicted EW variations 
from our "Toy" model and the observed EW variations is 
shown in Fig. [14] As illustrated, the model is able to quali- 
tatively reproduce the general features of the observed EW 
variations including the sharp emission minima, the smooth 
variation during maximum emission, and the double peaked 
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nature of the variability, suggesting that the Ha emission 
does form in an optically thick, flattened distribution of cir- 
cumstellar plasma with a disc-like structure. However, our 
model predicts a lower emission level between phases 0.0 to 
0.25 and 0.75 to 1.0 and a higher emission level between 
phases 0.25 to 0.4 and 0.6 to 0.75 than observed, and there 
is a visible offset between the phases of emission minima; our 
model predicts the minima to occur at phases 0.23 and 0.77, 
not 0.25 and 0.75 as observed. To investigate these discrep- 
ancies, we conducted another search allowing the disc axis 
inclination to vary. Our results show that a = 88 =b 1° (and 
a corresponding R[ n = 1.0 -R* and R ou t = 1.6 R*) provides 
a better overall fit to the observations, but is inconsistent 
with our derived value of f3. This disc orientation is better 
able to reproduce the emission level throughout the rotation 



cycle and correctly predicts the phases of minimum emission 
as illustrated in Fig. 1141 

In addition to investigating the magnetospheric prop- 
ertie s via the "To y " mo del, we also used the method 
of [S undqvist eT"al1 (l2Q12h to study the Ha variability. 
ISundqvist et al.l (|2012h showed that in a slow rotator such 
as HD 57682, the transient suspension of wind material 
in closed magnetic loops, which is constantly fed by the 
quasi-steady wind outflow, leads to a statistically overdense, 
low-velocity region in the vicinity of the magnetic equa- 
tor that causes persistent, periodic variations of Balmer 
li nes. This method utilises the 2D MHD wind simulations 
of £d -Doula. Owocki &; Towns end (2008) and the magnetic 
parameters determined from the Bn measurements to de- 
termine the properties of the magnetically confined plasma. 
We note however that the original MHD models were com- 
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Figure 8 - continued 



puted for an O-type supergiant. Therefore, in the Ha syn- 
thesis here, we have simply re-scaled the Ha scaling invari- 
ant Q — (M x fci) /(voo x R*) 3/2 (jPuls. Vink fe Naiarrol 
2008) from the original simulation to the value obtaine d 
for HD 57682 using the parameters of [Grunhut et all (2009). 
These properties are then used to synthesise the expected 
Ha line profile as viewed from our line-of-sight as the star 
rotates. Unlike the "Toy" model, these simulations are sen- 
sitive to the adopted mass-loss rate. We characterise the 
mass-loss rate as the product of the stellar mass-loss rate 
and the square-root of a wind clumping factor, which we as- 
sume to be about 4-10 as suggested by several recent studies 



.(em- 


Bouret et al.ll2003l:lNaiarro et alJl201ll:ISundavist et all 


boli 


). We attempted to constrain the wind properties by 



fitting the observed Ha profile during phases that have 
the lowes t contributions f rom t he disc emission. As already 
noted by I Grunhut et all (|20Q9h , we cannot fit the Ha pro- 
file using a mass-loss rate as derived from the UV lines 
(log M ~ -8.85 M© yr" 1 ). The additional emission resulting 
from the wind is too low at phases 0.25 or 0.75 (minimum 
emission) , which results in a line profile that is mainly in ab- 
sorption, contrary to what is actually observed (see Fig. ll5[) . 

Our results (based on fits to the Ha profile) suggest 
that the lack of emission is either due to an underesti- 
mated stellar mass- loss rate (by a factor of 10 to 30), or 
that the wind is very clumpy (note that clumping was n ot 
taken into account in the study of iGrunhut et al. I d2QQ§))- 
This is consistent with many recent studies that find large 
discrepancies between UV and optical mass- loss rates (see 
ISundqvist, Owocki &; Pulsl ([201 lh and references therein). 
Furthermore, according to Figs. [16] and [TTI we find that 
the model using the UV derived mass-loss rate produces 



little to no emission with any significant variability. In fact, 
the amplitude and the offset of the EW variations shown in 
Fig. [16] are very sensitive to the adopted mass-loss rate. If 
we choose a mass-loss rate that best fits the profiles at min- 
imum emission phases, we obtain a mass- loss rate of ~15 
times the suggested UV mass- loss rate (see Fig. I15|) . How- 
ever, while the predicted minimum emission of the EW vari- 
ations is well- fit, the maximum emission is underestimated. 
A mass-loss rate of ~20 times the previous estimate pro- 
vides a better fit to the EW variations, by underestimating 
the maximum emission and over-estimating the minimum 
emission. A mass-loss rate of ^30 times the previously es- 
timated value is needed to provide the correct amplitude of 
the EW variation, but, as the wind emission is predicted to 
be significantly higher, there exists a significant offset be- 
tween the predicted and observed EW values. In any case, 
our MHD simulati ons indicate t hat a higher mass-loss rate 
than derived by I Grunhut et al. (2009) is necessary in order 
to observe significant variability, with similar characteristics 
to the observed Ha dynamic spectrum (Fig. [TT)) . The char- 
acteristics of the MHD EW curves are similar to the predic- 
tions of our "Toy" model; the phasing of minimum emission 
are not at phases 0.25 and 0.75 and the relative amplitude of 
the emission is slightly lower when using a = 79° (note that 
the pseudo-disc naturally forms along the magnetic equator 
in these MHD simulations, i.e. a = j3). Just as was found 
with the "Toy" model, a higher f3 = 88° provides a bet- 
ter agreement with observations. Of particular interest, is 
the fact that the synthesised emission variation in Fig. [17] 
is predicted to be symmetric about the systemic velocity of 
HD 57682 with the adopted high f5. However, a lower f3 does 
provide increased occultation of the disc and therefore nec- 



© 2012 RAS, MNRAS 000.ITU2T1 



Investigating the spectroscopic, magnetic and circumstellar variability of HD 57682 13 



Hell 4686 



0.630 



g 0.640 
£ 0.650 




0.2 0.4 0.6 0.8 
Phase (P= 63.5708 d) 

Hel 5876 



0.700 
0.800 
0.900 

24.0 

"w 23.0 

B 22.0 

~" 21.0 
> 

x 20.0 
19.0 



l 1 1 1 l 1 1 1 l 1 1 1 l I 1 * l 1 1 1 l 1 1 1 l 1 1 1 



i. 



i i i i i i i i i i i i i i i i i i i i i i i i i i 



_L_ 



# • 
# 



§ 



# * 
i 



-0.2 0.0 0.2 0.4 0.6 0.8 
Phase (P= 63.5708 d) 



1.0 1.2 



2.300 



^ 2.400 
w 

2.500 



25.0 
20.0 

| 150 
^ 10.0 

25 5.0 

0.0 



I 1 1 1 I 1 1 1 I 1 1 1 I 4 1 1 I 1 1 1 I 1 1 1 I 1 1 1 



* i 



§ 



f I I I I I I I I I I I I I I I I I I * I f I I I I I I I 



-0.2 0.0 0.2 0.4 0.6 0.8 
Phase (P= 63.5708 d) 

Ha 



1.0 1.2 



G 

> 




-0.2 0.0 0.2 0.4 0.6 0.8 
Phase (P= 63.5708 d) 



1.0 1.2 



Figure 9. Rotationally phased EW measurements (top panels) and radial velocity measurements (bottom panels) for selected spectral 
absorption lines with clearly non-sinusoidal EW variations. Ha measurements from ESPaDOnS (black circles), CES (red squares), FEROS 
(blue Xs), Las Campanas Observatory (turquoise diamonds), UVES (pink up-facing triangles), and BeSS (green down-facing triangles) 
datasets. The BeSS Ha radial velocity measurements are shown as grey triangles. The large scatter of these measurements with respect 
to the other datasets likely results from their poor S/N and low spectral resolution. 



essarily results in a more asymmetric line profile variation 
at some phases, similar to the observations. 



7.2 Greater circumstellar region 

We continued our investigation of the circumstellar environ- 
ment to larger distances as probed at infrared wavelengths. 
Our original intention was to investigate the infrared prop- 
erties of HD 57682 's SED to address the uncertain classifi- 
cation of HD 57682 as a classical Oe star and to potentially 
probe the extent of the magnet osphere. 

In Fig. [HI we show the SED of HD 57682 from the ultra- 
violet to infrared wavelengths. The data were collected from 
the IUE database and the available photometric archives 
discussed in Sect. [2] Also included in this figure are the- 



oretical^ ED s based on the original stellar parameters of 
iGrunhut et al . (2009) and updated values presented here . 
A comparison between our cmfgen ([Hillier Mi ller 1998) 
model with the newly inferred parameters and the IUE data 
is also presented in the lower panel of this figure, illustrating 
the agreement between the observations and the model. We 
no te that we have still a dopted the lower UV mass- loss rate 
of IGrunhut et al . (2009) in this cmfgen model. Our new 
m odel has revised the d istance and reddening parameters 
of IGrunhut et al. | (|2009l . d ~ 1.3 kpc, E(B - V) = 0.07), 
which were only derived from the IUE data and do not pro- 
vide a good fit to the photometric measurements blue- ward 
of 50 000 A (5/im). Our revised parameters are d ~ 1.0 kpc 
and E(B - V) = 0.12. The SED also shows a large dis- 
agreement between the IRAS and WISE data points, but 
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Figure 10. Phased variations in selected spectral lines with sinusoidal equivalent width variations. Plotted is the difference between the 
observed profiles and the profile obtained on 24 December 2010 (dashed-red, top panel), which occurred at phase ~0.75, corresponding 
to when the magnetic equator crosses the line-of-sight. 
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Figure 11. Same as Fig. [10] but for selected spectral lines with clearly non-sinusoidal EW variations or EW variations that appear 
anti-phased with the other lines. Plotted is the difference between the observed profiles and the profile obtained on 24 December 2010 
(dashed-red, top panel), which occurred at phase ~0.75, corresponding to when the magnetic equator crosses the line-of-sight. 



this reflects the variation in the aperture sizes of the dif- 
ferent instruments; the IRAS aperture is about five times 
larger and therefore includes more contribution of circum- 
stellar emission (as further discussed below), while the WISE 
measurements represent more of the true "stellar" flux, and 
appears consistent with the other infrared measurements. 
The observed SED does show a large disagreement with the 
theoretical SED for measurements red- ward of 1 x 10 5 A (10 
/im). The observed infrared emission appears to be increas- 
ing with increasing wavelength and not decreasing as pre- 
dicted by contributions from pure stellar flux, and therefore 
represents an additional source of emission. 

In order to better understand the characteristics 
of the infrared ex c ess w e utilised the procedure of 
llgnace &; Churcriwe 11 (2004) to model the free-free emission 
that would be emitted from a spherical cloud; Infrared excess 
is common amongst Be/Oe stars and is caused by free- free 
emission from the circumstellar disc. Using this procedure 



we find that we can qualitatively match the characteristics 
of the infrared excess by varying the properties of the cloud. 
A brief investigation of this model shows that the location of 
the infrared bump is mainly controlled by the temperature of 
the cloud (we find T ~ 20000 K to be sufficient), but the size 
of the bump is controlled by the relative sizes of the cloud 
and star. For simplicity we assume that the cloud is of con- 
stant density, which then requires a cloud of approximately 
100 R* to reproduce the observations. If we instead adopt 
a cloud with many large clump s of different densities and 
at different radii (as is done in llgnace fe Churchw ell 2004) 
we can achieve similar results without the need for a 100 R* 
cloud. However, if we assume that the cloud is not highly 
clumped and results from stellar mass lost by HD 57682, this 
would imply an unreasonable mass for the cloud, suggesting 
that HD 57682 is in fact not a classical Oe star. 

In light of this conclusion we further investigated the 
IRAS photometry. In the end we suspect that the additional 
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Figure 12. Same as Fig. 1101 but for spectral lines that show evidence of two emission features per rotation cycle. The Balmer lines 
are plotted as the difference between the observed profiles and a NLTE tlusty model profile (top panel, dashed-red) to highlight the 
circumstellar emission. The He I A5876 lines is plotted as the difference between the observed profiles and the profile obtained on 24 
December 2010 (dashed-red, top panel). 
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Figure 13. Illustration of our "Toy" model for the Ha emission disc. The left panel provides an example schematic diagram showing 
the orientation of the plane of the disc for a given disc inclination a relative to the rotation axis. The other panels represent projections 
of the disc (solid black) and central star (solid red) onto our line-of-sight during phases 0.0, 0.33, and 0.66. The thin green arrow in these 
panels represents the projected rotation axis relative to our line-of-sight (taken to be 60°), while the thick blue arrow represents the disc 
axis of 79° relative to our line-of-sight. 
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Figure 14. "Toy" model compared to phased Ha equivalent 
width measurements (black circles). The solid red curve repre- 
sents a model with an a = 79° between the rotation axis and the 
disc axis, as inferred from the fits to the B% variations, while the 
dashed blue curve represents the best-fit model with an cm = 88° . 
Both models assume i = 60°. 
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Figure 15. Comparison between observed Ha profile (solid 
black) from 2010 December 24, corresponding to a phase of 0.72, 
and sy nthetic profiles following the procedure of lSundqvist et all 
(2012) at this same phase. The model profiles were created with 
an adopted mass-loss rate as determined bv lGrunhut et al.l (2009) 
(10 -8-85 ; dotted blue), a mass-loss rate of 17 (solid green) and 21 
times this value (dashed red.) 



emissio n is a result of H D 57682 's illumination of the IC 2177 
nebula (|Halbedellll993h . A look at the HIRES/IRAS maps, 
as shown in Fig. 1191 shows significant infrared emission cen- 
tred on HD 57682 at 25 /xm and 60 /im, with the size of the 
emission cloud larger at 60 /im. Emission is also evident at 
12 /im, but does not reach the same spatial extent as found 
at 25 or 60 /xm. At 100 /xm the cloud overwhelms the field of 
view (30 x 30 arcmin) , but appears to be a blend of several 
clouds. It is also evident from these images that the major- 
ity of the point sources in the field of view are illuminating 
the nebula at these wavelengths, further corroborating our 
suggestion that the infrared emission does not result from 
a circumstellar disc or even the extended magnetosphere of 
HD 57682. 
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Figure 16. Comparison between the measured EW variations 
of HD 57682 (black circles) and the theoretical EW variations 
of our MHD model for different mass-loss rates as labelled (UV 
corres ponds to the UV mass-loss rate derived by Gru nhut et al.l 
(2009)). While the absolute EW variations of the M/M = 30x 
UV (dash-dotted blue) model does not fit the observed EW curve, 
the amplitude of the predicted variations from this model does 
provide a good fit, as illustrated for the same model but vertically 
shifted to better match the variations implied by the observations 
(black dotted line). 



8 DISCUSSION & CONCLUSIONS 

In this paper we report on our ongoing work towards un- 
derstanding the magnetic properties and variability of the 
magnetic 09 subgiant star HD 57682. Our current analy- 
sis is based on an extensive dataset spanning almost 15 
years, including archival spectra from ESO's CES, UVES 
and FEROS spectrographs, and from the echelle spectro- 
graph at LCO. In addition numerous low-resolution spectra 
were acquired and accessed from the BeSS database. Fur- 
thermore, we also present 13 newly-obtained, high-resolution 
spectropolarimetric observations acquired with ESPaDOnS 
at the CFHT within the context of the MiMeS project. 

A period analysis performed on the spectroscopic Ha 
EW variations and the polarimetric Bi measurements re- 
sulted in a single, consistent period between the polarimetric 
and Ha variations of 63.5708 ±0.0057 d, which we infer to be 
the rotational period of this star. T his long rotational p eriod 
suggests that the v sin i inferred by I Grunhut et al.l (|2009h is 
too high by a factor of ~2, which led us to attempt to remea- 
sure this value. We found that the C IV A5801 photospheric 
line was the least susceptible to additional macroturbulence 
in the line profile and therefore were able to estimate a new 
vsini — 4.6 =b 0.6kms _1 , which is now consistent with the 
~63d period. This period also implies that the inclination 
of the rotation axis i > 30° with a best-fit value i ~ 56^26 °- 

From fits to the rotationally phased Bi measurements 
we were able to constrain the magnetic dipole parameters of 
HD 57682, assuming that the magnetic field is well char- 
acterised by the ORM. Taking i = 60°, the data imply 
that B d = 880 ± 50 G and f3 = 79 ± 4°. However, direct 
modelling of the individual mean LSD profiles suggests a 
somewhat weaker dipole field strength of ^700-900 G and 
a smaller obliquity angle of ^60-68°. However, the results 
are strongly dependent on the adopted depth and width of 
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Figure 17. Rotationally phased variations in obse rved (right pan e l) and modelled Ha profiles. The model in the left panel shows the 
predicted variations assuming the mass-loss rate of Grunhut et al . (2009), while the central pa nel uses an ado pted mass-loss rate of 21 
times this value. Note the lack of predicted variations of the model adopting the mass-loss rate of lGrunhut et al.l(|20Q9h and the predicted 
emission symmetry. 



the modelled LSD Stokes I and V profiles, as discussed in 
Sect. [5] Therefore, it is likely that improved modelling of 
the observed line profiles that correctly takes into account 
the profile variations could reduce the disagreement and we 
therefore suggest that the real magnetic dipole parameters 
are closer to the parameters implied by the Bg variation. 

The profile variations of nearly all observable lines in 
the high-resolution ESPaDonS spectra were also investi- 
gated (Figs. [8] and [9}. Our analysis suggests that we can 
essentially separate all lines into one of two categories as 
determined by the rotationally-phased EW variations of the 
spectral line. There are those lines that show clear evidence 
of double-peaked variations similar to Ha, while the rest 
show single-peaked, sinusoidal variations. The LPV of these 
lines, as described in Sect. [6] show rotationally phased veloc- 
ity variations of the pseudo-absorption or pseudo-emission 
components, and the pattern of this variability is generally 
consistent amongst all the lines. 

In Sect. we clearly showed that the double-peaked 
variations of the Ha EW is well explained by the presence 
of the variable projection of a flattened distribution of mag- 
netospheric plasma trapped in closed loops near the mag- 
netic equatorial plane. This flattened distribution is also 
predicted by MHD simulations using the st ellar and mag- 
netic properties as determined here and by I Grunhut et al.1 
(2009), although we find that a mass-loss r ate that is 20- 
30 tim es greater than the value inferred bv I Grunhut et all 
(2009) from the UV lines is necessary to produce the level 
of emission that is observed. However, at present we are not 
able to model UV resonance lines using our MHD wind sim- 
ulation and therefore must resort to modelling these lines 
using a spherically symmetric wind model. The predicted 
UV spectrum, as modelled using cmfgen, indicates the pres- 
ence of several intense spectral emission lines that are ob- 
served mainly in absorption in the IUE spectrum (e.g. Si IV 
A1400 and Niv A1718) or that the observed lines are too 
weak (e.g. Nv A1240 and Civ A1550). However, the IUE 
spectrum is well reproduced when adopting the lower mass- 



loss rate, as already demonstrated bv I Gru nhut et al.1 (|20Q9h 
and illustrated in the lower panel of Fig. 1181 Note that cm- 
fgen assumes a spherically symmetric wind, whereas the 
MHD wind simulations predict a highly non-spherical wind 
with important consequences for spectral line d iagnostics 
(jud-Doula fc Qwockil[200i : ISundavist et alJl2012h . Further- 
more, the UV lines are very sensitive to the ionisation struc- 
ture of the wind, which may also be largely affected by a 
non-spherically symmetric wind. This is not the case for Ha 
as it is a recombination-based optical emission line, which is 
insen sitive to the exact mod elling of the ionisation structure 
(e.g. [Sundqvist et al.l [2012). Therefore it is not surprising 
that there is a disagreement between the mass-loss rates in- 
ferred by our optical and UV analyses. Future MHD UV line 
modelling will be necessary to resolve this discrepancy. 

The MHD model provides not only an independent 
analysis of the stellar mass-loss rate, but also of the magnetic 
geometry, which is constrained by the characteristics of the 
Ha EW variations. Using the higher mass-loss rate shows 
that the Ha emission variation can be reasonably modelled, 
however there are several inconsistencies between the models 
and observed data that reflect our uncertainty in the mass- 
loss rate and potentially the magnetic geometry. Our results 
show that the inclination of the disc relative to the rotation 
axis must be greater than the inferred magnetic obliquity 
(as measured from our polarimetric data), since a higher 
disc inclination provides a much better fit to the observed 
EW variations. This may also be attributed to warping of 
the disc, or may even reflect the fact that the disc does not 
necessarily lie in the magnetic equatorial plane as predicted 
by MHD simulations. 

As with the LPV for the sinusoidally varying lines, there 
also appears to be rotationally phased velocity variations in 
the position of the central emission component in the Ha 
line profile that is not predicted by the MHD simulations. 
We discuss a few potential explanations. If the disc struc- 
ture is not vertically symmetric along the rotation axis (i.e. 
not symmetric about the equatorial rotation plane) the re- 
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Figure 18. Upper: Observed IUE observations and infrared 
flux measurements from the indicated sources. Also included 
are stella r models co rresponding to the parameters derived by 
iGrunhut etHI (120091 ) (d = 1.3 kpc and E(B - V) = 0.07; dashed 
grey) and the revised model presented here (d = 1.0 kpc and 
E(B - V) = 0.12; dashed red). The larger flux of the IRAS data 
points compared to the WISE measurements reflects the signifi- 
cantly larger aperture and therefore increased contribution from 
the surrounding emission. Lower: Comparison between the ob- 
served IUE spectrum (solid black) and the CMFGEN model with 
the parameters inferred in this work (dashed red) in the wave- 
length range of 1380-1 730 A. The cmfgen m odel still adopts the 
lower mass-loss rate of IGrunhut et all (2009) and not the higher 
value needed to fit Ha (Fig. [15] see Sect. [7] for further discus- 
sion). Both the model and observed spectrum were smoothed for 
display purposes. 



suiting line profiles would be asymmetric and would result 
in velocity shifting of the disc emission. Likewise, enhanced 
emission from a non-uniform velocity flow where the veloc- 
ity flow occurs at high magnetic latitudes could also ex- 
plain our observations. For example, during phases where 
the disc is viewed nearly face on, the emission could ap- 
pear systematically offset from the systemic velocity if there 
is an additional outflow or inflow, for instance, due to the 
stellar wind or circumstellar plasma falling back onto the 
star. This would result in no velocity offset during phases 
where the disc is viewed nearly edge on as the additional 
outflow or inflow would appear perpendicular to our line-of- 
sight. However, it is difficult to understand a scenario where 
there would be a preference for outflow or inflow of plasma 
at one magnetic hemisphere compared to the other. This 
is required to explain the velocity shifting of the emission 
that occurs as our line-of-sight rotates from one magnetic 
hemisphere to the other. Another possible solution is to in- 
troduce an offset of the dipole relative to the centre of the 
star. Our preliminary tests suggest that a small offset of 




Figure 19. Infrared emission as detected by IRAS/HIRES at 
different wavelengths (top left - 12 /im, top right - 25 /im, bottom 
left - 60 /im, bottom right - 100 /im). Each image represents a 
30 x 30 arcminute scale. Point sources are indicated by the red 
circles and the field of view is centred on HD 57682. Note that 
the elliptical shape of the point sources is a result of the scanning 
direction of the IRAS detector. 



~0.2i?* along the rotation axis could explain the velocity 
shifting of the emission that is observed in Ha, although 
MHD models implementing this offset are necessary to con- 
firm this phenomenon, which will be the subject of future 
studies. 

In addition to the LPV, we also measure RV variabil- 
ity in nearly all spectral lines and find that all lines with 
sinusoidally varying EW measurements show approximately 
sinusoidally varying RV measurements. iTurner et al.l ([2008) 
investigated the binarity of this star and found no evidence 
of a companion in either their /-band adaptive optics or RV 
measurements. If we assume that the RV variations are due 
to a binary companion, a least-squares analysis of the RV 
measurements obtained from the mean LSD / profile would 
suggest a mass function /(M/M©) = 2.2 ± 0.9 x 10" 6 for 
the companion star. For reasonable orbital inclinations, this 
would require the companion to be a very low-mass (and 
presumably a low-luminosity) star, and therefore should not 
have a significant effect on the observed line profile, which we 
clearly do observe. Taking this into account, it is therefore 
unlikely that the RV variations are due to a binary compan- 
ion since we do not observe clear shifts of the entire spec- 
troscopic lines. We propose that the RV variations are likely 
a result of the variable asymmetry of the line profiles and 
simply reflects the varying location of the centre-of-gravity 
in the profiles. This is clearly evident in the lines that show 
a double-wave EW curve since the RV measurements also 
show a more complex double- wave pattern. Furthermore, 
there also appears to be more of a direct relationship be- 
tween between the RV and EW measurements in these lines, 
while the majority of the sinusoidally varying lines show a 
^0.25 phase offset between the RV and EW variations. 

The root cause of the LPV in the sinusoidally vary- 
ing lines is still an outstanding issue. It is highly unlikely 
that these variations are caused by surface features (e.g. 
chemical "spots") as the observed variability between mul- 
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tiple ions of a given species does not behave as predicted 
by LTE or NLTE spectrum synthesis. For example, our ob- 
servations demonstrate that the weak Niv A4058 line ex- 
hibits the same relative EW variation as the strong Nlll 
A4523 line, but spectral modelling of these lines show that 
the weaker line should show 30 percent stronger variation 
if due to abundance changes. Furthermore, the strong Si IV 
A4115 line and the weak Si IV A4403 line show nearly the 
same level of EW modulation, but the weak Si IV line should 
show a 60 percent larger amplitude in the EW variation if 
resulting from abundance changes. The weak Si IV A4667 
line is also predicted to show the same relative EW varia- 
tion as the A4403 line, but our observations show a signifi- 
cantly greater level of variability (about 50 percent larger) 
in the A4667 line. Moreover, abundance spots require that 
all spectral lines of a given species vary in phase. This is 
clearly not the case for the Civ lines (5801 A, 5811 A) that 
vary in antiphase with the C in lines (as shown in Fig. [8]), or 
the Niv A4058 and N,iii A4523 lines, which appear shifted 
relative to one another by a quarter of a cycle. The inter- 
pretation of the variability in terms of surface spots would 
also imply that all elements (except perhaps C) are con- 
centrated at the same general region on the stellar surface, 
which would be a remarkable new phenomenon. Another 
possibility is that the variations may result from pulsations; 
however, pul sations should occur on much shorte r time - 
scales (e.g. lAerts, Christensen-Dalsgaard &; Kurt j l2010h . 
The rotationally-phased, double-peaked EW variations are 
well explained by the variable projection of a flattened dis- 
tribution of magnetospheric plasma, and at first glance it is 
difficult to apply this same disc model to explain the vari- 
ations in the other lines. Occultation or variable emission 
from a symmetric disc would result in double-peaked vari- 
ations in all lines and not the single-peaked emission that 
we observe. However, the LPV of the higher Balmer lines 
(Fig. I12[) likely provide the link between the double-peaked 
variability and single-peaked variability; a decrease in the 
strength of the emission feature at phase 0.25 and 0.75 com- 
pared to the emission feature at phase is clearly observed 
in higher Balmer lines. Furthermore, this is also illustrated 
in the EW variations of the Balmer lines (Figs. [8] and [9J, 
which show a clear transition from the double-peaked EW 
curve at Ha to a single-peaked curve at H7. Therefore, we 
tentatively conclude that all line profile variability observed 
in HD 57682 is a result of variations of the flattened distri- 
bution of magnetospheric plasma. The details of the vari- 
able emission for different lines likely reflects the fact that 
each line probes a different region of what is likely an asym- 
metric spatial or velocity distribution of the magnetospheric 
plasma and that the disc is not completely optically thick 
for all lines-of-sight, for all wavelengths. We suspect that full 
3D MHD models and an extension of th e line profile synthe- 
sis techniques of ISundqvist et aL I (|2012h to lines other than 
Ha are necessary to properly model these properties. 

In summary, the narrow-lined photospheric spectrum 
and weak wind of HD 57682 have allowed us to study its 
magnetic and magnetospheric properties at a level of de- 
tail not currently achievable for any other magnetic O-type 
star. The results of our analysis indicate a highly com- 
plex behaviour that can be used as a testbed for future 
more sophisticated 3D MHD simulations to better under- 
stand non-rotationally supported magnetospheres (also re- 



fer redtoas dynamical magnetospheres IPetit et al.1 l201ll : 
ISundqvist et al.ll2012r ) of hot stars. 
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Table Al. Spectroscopic Ha EW measurements. Included is the 
instrument or observatory name, the heliocentric Julian date, and 
the measured Ha equivalent width (EW) and its corresponding 
la uncertainty. 



Instrument 


HJD 


ria. EW (A) 


cr EW (A) 


CES 
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0.007 


CES 
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0.008 
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0.735 


0.012 


CES 
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0.569 


0.011 


CES 


2450130.6194 


0.454 


0.010 


CES 


2450131.6155 


0.311 


0.007 


UVES 


2452621.8553 


0.278 


0.035 
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2453042.6286 


1.130 


0.033 


FEROS 


2453043.6016 


1.468 


0.029 


FEROS 


2453045.6764 


2.187 


0.066 


FEROS 


2453045.6873 


2.154 


0.034 


FEROS 


2453361.6916 


1.308 


0.040 


LCO 


2455111.8964 


1.612 


0.015 


LCO 


2455112.8549 


1.701 


0.019 


LCO 


2455113.8919 


1.560 


0.019 


BeSS 


2454786.6621 
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0.032 
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1.517 


0.074 
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2455233.4676 


1.259 
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1.341 


0.045 
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1.667 


0.083 
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2455235.4449 


1.630 


0.094 


BeSS 
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2.019 


0.041 
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0.065 
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2455240.4484 


1.816 


0.058 
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1.908 


0.082 
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1.205 


0.059 
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0.140 
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0.121 
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0.137 
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0.103 
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0.061 


BeSS 
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0.075 


BeSS 
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0.102 
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0.301 


0.114 
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